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Abstract 

This paper presents an extended Kalman filter for estimating attitude sensor timing errors. Spacecraft attitude is determined 
by finding the mean rotation from a set of reference vectors in inertial space to the corresponding observed vectors in the body 
frame. Any timing errors in the observations can lead to attitude errors if either the spacecraft is rotating or the reference 
vectors themselves vary with time. 

The state vector here consists of the attitude quaternion, timetag biases, and, optionally, gyro drift rate biases. The filter 
models the timetags as random walk processes: their expectation values propagate as constants and white noise contributes 
to their covariance. Thus, this filter is applicable to cases where the true timing errors are constant or slowly varying. 

The observability of the state vector is studied first through an examination of the algebraic observability condition and then 
through several examples with simulated star tracker timing errors. The examples use both simulated and actual flight data 
from the Extreme Ultraviolet Explorer (EUVE). The flight data come from times when EUVE had a constant rotation rate, while 
the simulated data feature large angle attitude maneuvers. The tests include cases with timetag errors on one or two sensors, 
both constant and time-varying, and with and without gyro bias errors. 

Due to EUVE’s sensor geometry, the observability of the state vector is severely limited when the spacecraft rotation rate is 
constant. In the absence of attitude maneuvers, the state elements are highly correlated, and the state estimate is unreliable. 
The estimates are particularly sensitive to filter mistuning in this case. The EUVE geometry, though, is a degenerate case 
having coplanar sensors and rotation vector. Observability is much improved and the filter performs well when the rate is 
either varying or noncoplanar with the sensors, as during a slew. Even with bad geometry and constant rates, if gyro biases 
are independently known, the timetag error for a single sensor can be accurately estimated as long as its boresight is not too 
close to the spacecraft rotation axis. 


1. Introduction 

One important step in spacecraft telemetry processing is assigning timetags to each sensor measurement. Errors in 
these timetags can occur as constant biases due to inaccurate time offsets used in the software or as possibly time-vary ing 
errors due to problems with the sensor itself or with input/output to an onboard data buffer. Timetag errors can lead to 
errors in the estimated attitude if either the spacecraft is rotating or the reference vectors themselves vary with time. 

This work examines the feasibility of estimating attitude sensor timetag errors along with the attitude quaternion and 
gyro rate bias (other calibration parameters could be included as well). The timing errors are thought of as constant biases 
(offsets from the true observation time), but the method also applies to slowly varying errors. 

The timetag biases are appended to the system state vector and estimated using an extended Kalman filter. Section 2 
describes the sensitivity matrix and state transition matrix needed as input to this algorithm. Some insight into the obser- 
vability of the larger state vector is obtained by examining the observability matrix in Section 3. Results of several test 
cases are presented in Section 4. The tests include examples with timing errors on either one or two sensors, with and 
without gyro rate biases, for timespans with and without attitude maneuvers. Examples where the true bias is time-varying 
are also considered. Conclusions are given in Section 5. 


* This work was supported by the National Aeronautics and Space Administration (NASA)/Goddard Space Flight Center 
(GSFC), Greenbelt, Maryland, under Contract NAS 5-31500. 


379 



2. Theory 


The extended Kalman filter theory, as applied to attitude and gyro bias estimation, is given in Reference 1. One dis- 
tinguishes the full state vector consisting of the quaternion and gyro bias, from the error state consisting of a small angle 
rotation vector, a, and a correction to the gyro bias, A b. For each sensor measurement, the resulting error state updates 
the prior estimate of the full state: it rotates the attitude and additively corrects the gyro bias. Between measurements, the 
full state is propagated using bias-corrected gyro data. 

This approach moves the normalization constraint on the attitude quaternion onto the 4 th component of the error 
quaternion. (The error quaternion is the quaternion corresponding to the rotation vector a . It automatically has 4 th com- 
ponent equal to unity, to first order in a.) This leaves the three components of a completely unconstrained, and the 
normal extended Kalman filter theory then is applicable. 

References 2 and 3 describe the unit vector filter (UVF) observation model used here. In this model, the actual 
measured quantities are converted into unit vectors in the body frame before being handed to the Kalman filter. The 
effective sensor noise is taken to be isotropic; that is, a single scalar parameter represents the strength of the noise, 
regardless of the direction of the unit vector in the actual sensor field of view. 

The rest of this section expands the UVF to include timetag biases by developing the sensitivity matrix and the state 
transition matrix. The sensitivity matrix relates the unit vector observations to the error state. It is needed at each sensor 
measurement to construct the Kalman gain and estimate a new error state. The state transition matrix is used to propagate 
the covariance of the error state between observations. 


Sensitivity Matrix 


The observation of a vector quantity (eg., a star unit vector) is modeled as a reference vector, v, rotated into the body 
frame at the true time of measurement, plus random, zero-mean sensor noise, n k . The attitude matrix, A true> rotates vectors 
from the inertial frame to the body frame representation. The predicted observation uses the attitude estimate conditioned 
on the prior observations 


true observation = w. = /f v, + n. 

true true true k 

predicted observation = w^ k _ { = A k \ k _ x v k 


( 1 ) 


The k\j notation indicates a quantity estimated using observation data through time t jy and propagated (if k*j) to time t k . 
Thus, in the usual parlance, k\k-\ indicates an a priori estimate and /r|A: an a posteriori estimate. 

The observation is assigned timetag t k , whereas the true time of the measurement is t^ = t k - x. In particular, t is 
the time bias to be estimated. 

The true timetag bias, its estimate, and the time bias error are related as x = + A x . At each update, one estimates A x 

and adds it to x^ to get the a posteriori x^. This propagates as a constant to the next update time, becoming x k ^ k . 

To simplify the notation, one can absorb the prior estimate x^j t k . The reference vectors always are obtained 
for corrected times using the best current time bias estimate, so the subscript-/: properly should indicate variables at that 
corrected time. Specifically, t k is corrected and replaced hereafter by t k -z^ k l9 and the time bias error becomes 

At ( 2 ) 
The advantage is that only At (and not need appear in the following derivation. 
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( 3 ) 


The attitude at (corrected) time t k is related to the attitude at the true observation time by 




where R A ( At) is the rotation matrix for any spacecraft motion during time At. One can represent this rotation by 




e 


-ia A Atx] 


(4) 


where the true rate, is the gyro output, u k , corrected for the gyro bias 

$ k = u k - b k (5) 


The estimate (S^j = u k ~ b k ^ k _ { is used to approximate ti k in Equation (4). Also, for any vector v, the ’’cross-product 
matrix" is defined 


0 


[vx]s 



( 6 ) 


The reference vector also may vary during the time interval. Its magnitude does not contribute information to the 
state estimate, so only rotations need be considered. Thus, the reference vector at time t k is related to the reference at the 
true observation time by 

v; = R v (Ax) v true (7) 


If the reference vector is rotating at a mean rate r k , then 


* v (At) = 


IF^AtxJ 


( 8 ) 


The sign in the exponent here is opposite to that in Equation (4) because R v refers to actual rotation of a vector rather than 
its apparent motion due to rotation of the body frame. 


The sign of the attitude error angle a is defined such that the true and estimated attitude matrices at time t k are 
related by 


A 


k 



(9) 


The seven-component state error vector is 



At 


k 


( 10 ) 
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The sensitivity matrix, H k , is the derivative of the observation with respect to the state, evaluated at the prior state 
estimate. One can also write this as a relationship between the observation residual, y k , and the error state 

y* = H k x k + < n > 


Note that the a priori error state is identically zero since it already has been used to update the full state vector after the 
previous observation (the actual error may be nonzero, but its estimate vanishes). Expanding y k to lowest order in the error 
components (expanding about zero) yields 

ft = - *Vi 

= A . v,_ + n . - A 


true’ true tw k **k\k- l^Jfc 

" (/ + [w t AT t x])(/ + ta Jt x])^ 1|;k . 1 (/-[r t AT t x])v t - + n k 

“ -Ki*-i*i s * - ~ + 

= *] At* + «* 


(12) 


Thus, the sensitivity matrix is 

H k = X 1 ^3x3 

If one simultaneously observes two vectors, vv 1 and w 2 , ( e.g using two star trackers) to estimate the timetag bias of only 
w 1 , one obtains 




-[wSc] 0 3x3 -vv ] x((5 -Af) 


- 2 - 1 0, 


- w x 


3x3 


(14) 


where all quantities are a priori estimates (k\k-l). 
State Transition Matrix 


The time evolution of the attitude and gyro bias parts of the error state is described in References 1 and 3. The 
attitude error is rotated at the rate o>, with a correction due to the gyro bias error and a noise term. The gyro and timetag 
biases are represented using random walk models. The timetag bias may be more properly modeled as a random constant, 
but the random walk leads to a somewhat more robust filter and prevents the covariance from underflowing. The evolution 
equation is 


dx r 
— = Fx + n 
dt 


[<3 x] / 3x3 0 3x/> 

a 



^3x3 ^3x3 ®3xp 

A b 

+ 


0 - 0,0 

pxi px 3 p*p_ 

Af 


V 


where Af here is written more generally as a p-component vector, representing timetag bias errors for p different sensors. 
The vector, r \ , is the white noise source driving random walks in the attitude, gyro bias, and timetag bias. 
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Assuming the angular rates are constant during the propagation time interval, At = t 2 - r p the transition matrix is the 
exponential of the evolution matrix F times At: 



4> 

* 

0 3 xp 

4>(f 2 ,f,) s 4>( AO = e™ = 

^3x3 

Ax3 

O 

U) 

4 


^px3 

^px3 

Ipxp. 


(16) 


with 


rn = ^(-[«x]y 

h o J- 


(17) 


and 


♦(0 = EH«x])' 


^1! 


(18) 


These sums collapse into well-known, finite expressions for <J> and i|r (eg., Reference 4). Finally, the transition matrix, 
along with the process noise ^-matrix obtained from r \ , is used to propagate the state error covariance over the time 
interval At (References 4 and 5). 


3. Observability 


The observability condition for an border, constant coefficient, deterministic, discrete system is that the matrix 


K = 


H 

//O 

//O 2 




(19) 


be of rank n (Reference 5). If there are two vector observations (six components) per update as in Equation (14), IK will 
have dimension 6nxn; the rank n requirement then is that there be an nXn nonsingular submatrix (/.e, its n remaining rows 
or columns must be independent after deleting 5 n of its 6n rows). 

This can be related to the batch solution for the epoch value of a state vector, x o . The state is observable if x o is 
uniquely determined by a set of n distinct observations, y k> made on x o and the propagated states 0 *jc o . From Equations 
(11) and (19) one has 


= <Kx ( 2 °) 

o 


>0 

tl 
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where the noise term has been dropped (it is not usually important for observability). Left multiplication by 0f r IF, where W 
is an arbitrary, positive-definite weight matrix, yields an nxn square matrix, yC^WK, on the right hand side. Now, if H 
has full rank, < K T W'K can be inverted to obtain 

*> 

x 0 = ( < H T W , H)~ l M r l V ( 21 ) 

?n- 1. 

and x 0 is recovered. 

When the attitude, gyro bias, and a timetag bias for a single sensor (p = 1) are estimated, there are n = 7 error state 
components, and one obtains from Equations (14), (16), and (19) 

[w 1 x] 0 w 1 xc3 

[w z x] 0 0 

[vv 1 x](J) [VP 1 x] ip VV 1 XG> 

[w 2 x]< J> [vv 2 x]i|j 0 

[w^x]^ 2 [>v 1 x](<J>x|r ^^1^) vv 1 xd) 

n-2 

[w^x]^" -1 [H? 1 ^ <^4* vv 1 xg5 

;=0 

n-2 

[vv 2 x](J) rt 1 [w 2 x]J}^t 0 

7=0 

where the reference vector time derivative, r , has been assumed negligible* 

Note that the assumption of a constant coefficient system was made to keep Equation (19) simple. Thus, the vectors w l 
and vv 2 were taken to be the same after each propagation. But if they are constant in the body frame, they must refer to 
different inertial frame vectors (assuming <}>*/). So, Equation (22) represents the state observability after measuring n 
distinct vector pairs. 

To make sense of Equation (22), evaluate it explicitly for a set of vectors with a geometry similar to that for the 
Extreme Ultraviolet Explorer (EUVE). Let 

d> = ( (*) , 0, 0 ) T 

w'=( 1,1,0 ) r /v/2 ( 23 > 

w 2 = (-1,1,0 flfi 
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and assume A t is small to simplify <J> and i|f Then, 


0 0 0 
0 0 0 
0 0 0 


0 

0 



0 1 

0 -1 

1 0 

0 1 

0 1 

-1 0 

-wAf 1 

u)A t -1 

1 w At 

-Cl) At 1 

-G) At 1 

-1 -coAf 

-2oA t 1 

2oAf -1 

1 2o>Af 


0 0 0 
0 0 0 
0 0 0 

0 0 At 

0 0 -At 

-At At 0 

0 0 At 

0 0 At 

-At -At 0 

0 0 2 At 

0 0 -2Af 

-2Af 2Af 0 


-G) 

0 

0 

0 

0 

0 

— G) 

0 

0 

0 

0 

0 

-G) 


(24) 


Several features can be identified. First, if one is solving only for the attitude {n = 3), only the leftmost 18x3 submatrix 
applies. One can see that the first, third, and sixth rows are independent and span the three-component vector space. Thus, 
as expected, a single pair of vectors is adequate to determine an attitude (it is not necessary to use rows beyond the sixth, 
so additional observations are redundant). 

Second, if one is solving for attitude and gyro bias (n = 6), the leftmost 36x6 submatrix applies. Again, rows 1, 3, 
and 6 span the attitude part of the state. Then rows 7, 9, and 12 span the gyro bias components. So this state is observable 
after two distinct observations of two vectors each. 

Third, if one is solving for attitude and time bias (« = 4), the 24x4 submatrix consisting of the leftmost three columns 
and the rightmost column applies. One can span the first three components, with no admixture of the time bias column, 
by using rows 1, 6, and 7. Then, row 3 completes the set. 

Finally, the entire 42x7 matrix has only rank 6. There do not exist seven independent rows. (This holds with the 
full cj) and , not just for the small At approximation.) The seven-component error state is not observable under the 
assumptions of a constant coefficient system and the EUVE-like geometry given in Equation (23). One can estimate the 
attitude and either the gyro bias or the timetag bias but not all three. 

In Equation (23), the rotation rate and observed vectors in the body frame are coplanar. The unobservability of the 
seven-component system actually arises from this degeneracy. If they are noncoplanar, ‘H does have the full rank 7. In 
addition, the state is observable if the system is time-dependent. Then, different rates effectively appear in subsequent 
rows, enabling all the components to be disentangled. Thus, for EUVE, slewing the spacecraft can improve observability 
both by introducing a time-dependence and by shifting the rotation vector out of the plane defined by the sensors. 
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4. Test Cases 


This section presents test cases showing when a timetag bias can and cannot be determined, and how it affects the 
estimation of the rest of the state vector. Most of these results are summarized in Table 1. For all cases listed in the table, 
except where indicated, the state vector consists of attitude, gyro bias, and timetag bias. 

The EUVE flight data cover a two-orbit timespan for December 16, 1992, when the spacecraft was in Survey Mode, 
rolling about the body x-axis at 3 rotations per orbit (rpo). Performance statistics reported below are obtained as mean 
values or root-mean-square (RMS) averages over the second orbit of the data set. The simulated EUVE data cover a 4000 
sec timespan with characteristics similar to the flight data. Mean and RMS values for these cases are obtained from the 
second half of the timespan. 

Two Timetag Biases 

The UVF software was modified to solve simultaneously for separate timetag biases for the two fixed-head star 
trackers (FHST1 and FHST2). Simulated biases were added to the true timetags by the driver routine. Inclusion or exclu- 
sion of gyro biases in the state vector did not change these results significantly. It was found that the algorithm could not 
determine the time biases. 

Figure 1 shows the estimated biases, their error (square root of the estimated variance), and the attitude error (root 
of the trace of the attitude part of the covariance matrix). The true (simulated) time biases were -0.8 and 0.0 sec. The 
estimates were -0.4 and 0.4 sec. The error is much larger than the estimated la error of 0.06 sec. Also, the attitude error 
is 40 arcsec, substantially larger than the 5-arcsec errors typical of the UVF in the absence of timetag biases. In addition, 
the actual deviation between the estimated attitude and the onboard computer (OBC) estimate is 250 arcsec (RMS over 
the second orbit). The corresponding deviation is only 12 arcsec when no timing errors are present (Reference 4). This 
attitude deviation comes from the 0.4-sec average timetag error times the 3 rpo roll rate, which leads to a large roll offset. 
The correlation coefficients obtained from the orbit-averaged covariance matrix are near unity between the two timetags 
and also between the timetags and the roll angle. It is interesting that in this, and several similar test cases, the algorithm 
correctly obtained the difference between the FHST1 and FHST2 timetag biases to within about 0.05 sec. (This did not 
hold for differently placed sensors.) 

Single Timetag Bias without Gyro Bias 

In all subsequent tests, the filter solves for the timetag bias for only one sensor at a time. In this example, FHST2 
is given a time bias of -0.3 sec. The true gyro biases are compensated using values obtained from the UVF smoother (Ref- 
erence 4) and are not estimated. Figure 2 shows the estimated FHST2 timetag bias and the timetag and attitude errors 
derived from the covariance matrix. The time bias converges in a few hundred seconds and settles on a mean value dif- 
fering from the true bias by only a few microseconds. Its estimated lo error is 0.005 sec (the actual standard deviation 
of the estimate is about 0.002 sec). The attitude estimates are very good (similar to solutions with no timing errors): 
estimated lo error is 4 arcsec, and RMS deviation from the OBC estimate is 12 arcsec. The largest correlation coefficient 
occurs between the timetag bias and the roll error, but its value is only 0.43. This is the only example for which the time 
bias estimation worked well without an attitude maneuver. 

Single Timetag Bias with Gyro Bias 

Figure 3 shows a case identical to that in Figure 2, except that gyro biases are also estimated. Comparison with 
Figure 2 indicates that the convergence is much slower, and, in fact, the time bias is moving in the wrong direction for 
much of the last few thousand seconds. The true simulated time bias is -0.3 sec, but the mean estimate (for the second 
orbit) is -0.25 sec. The estimated lo errors obtained from the covariance matrix are 0.02 sec, 12 arcsec, and 0.03 deg/hr 
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for the time bias, attitude, and gyro biases, respectively. These are significantly worse than the 5 arcsec and 0.007 deg/hr 
attitude and gyro bias errors obtained without timetag errors. The RMS deviation from the OBC estimate is 26 arcsec. 
In addition, the correlation coefficients indicate that the solution is unreliable, several being over 0.9. These results show 
that EUVE’s attitude, gyro bias, and timetag bias together are not observable without maneuvers. 

Having established that the filter can correct either the timetag or the gyro biases, but not both, iterative methods were 
attempted. Data with both a gyro bias and an FHST2 timetag bias were processed repeatedly, solving alternately for the 
gyro and the timetag bias. The results do not converge. Similarly, the iterative solution for two separate FHST timetag 
biases fails to converge. Thus, as expected, the observability problems that spoil the simultaneous estimation of two time 
biases or a timetag and gyro bias cannot be circumvented by estimating them separately and iteratively. 

Timetag Bias Estimation During Maneuvers 

A series of tests using simulated EUVE data were performed to show the influence of attitude maneuvers on the state 
observability (attitude, gyro bias, and timetag bias). The simulations include 0, 1, or 2 slews of roughly 30 deg at 
0.1 deg/sec. The first slew is about the body _y-axis. When there are two slews, the second is about the body 2 -axis. Both 
slews occur during the first 2000 sec. 


One advantage of simulation data is that the noise characteristics are known exactly, so these tests also offer an 
opportunity to look briefly at the effects of filter mistuning. The filter is precisely tuned when it agrees with the simulator 
on the strengths of the white noise sources. These are the process noise q in Equation (15) and the sensor noise h in 
Equation (1). In these tests, the simulated timing error is a constant (q t = 0) while the filter’s timetag noise parameter is 
set to a small nonzero value. All other noise parameters were tuned exactly. In this case, the actual solution deviations 
from the true values agree very well with the expected lo errors. When the spacecraft undergoes maneuvers, these errors 
decrease in parallel. The lo timetag error decreases from 0.024 to 0.009 to 0.006 sec for 0, 1, and 2 slews, while the 
actual mean deviations go from 0.021 to 0.012 to 0.004 sec (see Table 1). The attitude errors and RMS deviations from 
the truth model both decrease from 14 to 7 to 5 arcsec. Most importantly, the largest correlation coefficient decreases from 
near unity to 0.77 for the two slew case, indicating that the separate state components have become distinguishable. 


When the filter is mistuned, the lo errors and actual deviations no longer agree so closely. For this test, the FHST 
noise parameter was set to 0.005 deg rather than the true 0.01 deg used in the simulator. This causes the expected errors 
to decrease while the mistuning should cause the actual deviations to increase. With no maneuvers, this is exactly what 
is found. The timetag deviation increases to 0.032 sec and the attitude deviation increases to 18 arcsec. What is interesting 
is that with two slews the solution is much less sensitive to this mistuning. In this case, the actual errors revert to near 
their tuned values, 0.004 sec and 5 arcsec. (Some degradation is expected, but this can only be seen by averaging results 
from many runs using different random number sequences and initial conditions.) One concludes that state observability 
is an important part of filter robustness. 

Estimation of Variable Timetag Bias 

For the final tests, the filter is asked to solve for a single, but time-varying, timetag bias. A sawtooth function with 
amplitude -0.5 sec and a period of either 60 sec or 3000 sec is added to the true FHST1 timetag. Gyro biases are not esti- 
mated; the rates are corrected using the true gyro biases, determined separately. 

If the filter memory time is too long, it will not be able to follow a varying time bias; it will converge instead to an 
average value. The process noise needs to be large enough to produce an uncertainty roughly comparable to the sawtooth 
amplitude in a time equal to the period. The timetag process noise is characterized by the parameter a . , obtained from 
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( 25 ) 


the two-time expectation of the white noise source ri t appearing in Equation (15): 

^[•n x (r) ■n T (f / )] = 

A range of values for a \ from 0.01 to 0.001 sec 2 /sec leads to similar results. (All other test cases where the timetag biases 
were constant used o\ = 10~ 8 sec 2 / sec.) 

Figure 4 shows the timetag bias from a typical 400-sec timespan for the more rapidly varying case. The estimate is 
shown as a solid line and the true timing error as a sawtooth dashed line. When the observability is good (stars available 
in both trackers), the estimated bias can follow the sawtooth with a residual of roughly 0.05 sec or less. However, the 
overall performance is not very good since the sawtooth quickly moves away from the estimated time bias between star 
identifications. 

In spite of the rather poor timetag estimation, inclusion of the time bias in the state vector does improve the attitude 
determination significantly. The RMS deviation from the OBC attitude is only 11 arcsec, as good as tests with no timing 
errors. Without time bias estimation, the sawtooth timing error causes an RMS attitude deviation of 74 arcsec. 

When the period of the sawtooth is increased to 3000 sec, the filter has much less trouble following it. Figure 5 
shows the timetag bias estimate (solid) and the slowly varying true bias (dashed). The RMS deviation from the OBC is 
11.5 arcsec in this case. 

5. Conclusions 

It has been shown that, without maneuvers, attitude sensor timing errors onboard EUVE and spacecraft with similar 
geometry can be solved for only for very limited circumstances. Examination of the Tf -matrix shows that the attitude, gyro 
bias, and timetag bias are not simultaneously observable. A Kalman filter solution leads to results with large correlation 
coefficients. These solutions are not robust: filter mistuning changes the results and iterative methods lead to radically 
different estimates or diverge altogether. The timetag bias can be accurately determined only if the gyro bias is not solved 
for and is compensated using an accurate independent estimate. 

Similarly, when there are timetag errors on two sensors, the filter algorithm cannot distinguish them. Large correla- 
tions between the individual sensors’ timetag biases indicate that their estimates are unreliable. Additional tests are needed 
to examine cases with three or more sensors, with and without attitude maneuvers. 

The observability improves dramatically when the spacecraft undergoes attitude maneuvers. The correlation coeffi- 
cients decrease and the state estimates improve for each slew. This only holds, though, with fairly large maneuvers. Tests 
with recent EUVE flight data having a pair of 2 deg slews showed no noticeable increase in observability. 


Even when the gyro bias has been compensated, the filter does not do well estimating a short-period, time-varying 
timetag bias. However, the problem here is probably insufficient data relative to the bias time scale (period equal to I 
minute in this test case). The filter has no trouble following the long-period variations in the final test. In either case, 
estimating the timetag bias improved the attitude estimates significantly. If a particular sensor has a known short-period 
timing problem, a better sensor model could be developed to estimate, for example, the period, amplitude, and phase as 
three new state vector elements, rather than trying to follow the time-variation using a single timetag bias with a large 
process noise. 
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Table 1 . Timetag Bias Filter Statistics 


Test Cases 

Mean Timetag 
Deviation 
from True 
(sec) 

Estimated 
Timetag 
la Error 
(sec) 

RMS Attitude 
Deviation 
from True 
(arcsec) 

Estimated 
Attitude 
1 a Error 
(arcsec) 

Largest 

Correlation 

Coefficient 

No Timing Error 

NA 

NA 

11.8 

5.0 

0.62 

One Timetag 
(w/o Gyro Bias) 

< 5 . E -06 

0.005 

12.0 

4.2 

0.43 

One Timetag 

0.052 

0.020 

25.8 

12.4 

0.97 

Two Timetags 

0.4 

0.06 

247 

40 

0.99 

Simulation with 
No Maneuvers 

0.021 

0.024 

13.5 

13.9 

0.98 

Simulation with 
One Maneuver 

0.012 

0.009 

6.9 

6.5 

0.85 

Simulation with 
Two Maneuvers 

0.004 

0.006 

5.1 

5.1 

0.77 

Simulation with 
No Maneuver • 
Mistuned Filter 

0.032 

0.010 

18.4 

7.4 

0.96 

Simulation with 
Two Maneuvers - 
Mistuned Filter 

0.004 

0.004 

4.8 

3.6 

0.65 
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Figure 1. Timetag Biases, Timetag Errors, and Attitude Error 
When Estimating Timetag Biases for Two Sensors 
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Figure 2. Timetag Bias, Timetag Error, and Attitude Error 
When Estimating Timetag Bias for a Single Sensor 
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Figure 3. Timetag Bias, Timetag Error, Attitude Error, and Gyro Bias Error 
When Estimating Gyro Biases and Timetag Bias for a Single Sensor 
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Figure 4. Estimated and True Timetag Biases When Estimating a 
Rapidly Varying Timetag for a Single Sensor 



Figure 5. Estimated and True Timetag Biases When Estimating a 
Slowly Varying Timetag for a Single Sensor 
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